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Abstract: 

A standard approach to analyzing tunneling processes in various physical contexts is to use instanton 
or imaginary time path techniques. For systems in which the tunneling takes place in a time dependent 
setting, the standard methods are often applicable only in special cases, e.g. due to some additional 
symmetries. We consider a collection of time dependent tunneling problems to which the standard 
methods cannot be applied directly, and present an algorithm, based on the WKB approximation 
combined with complex time path methods, which can be used to calculate the relevant tunneling 
probabilities. This collection of problems contains, among others, the spontaneous nucleation of 
topological defects in an expanding universe, the production of charged particle - antiparticle pairs in 
a time dependent electric field, and false vacuum decay in field theory from a coherently oscillating 
initial state. To demonstrate the method, we present detailed calculations of the time dependent decay 
rates for the last two examples. 
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1 Introduction 



Problems involving quantum mechanical tunneling in a time dependent setting can arise in 
a wide variety of contexts, such as the ionization of atoms by strong laser fields [jl|, pair 
creation of charged particles in time dependent background electromagnetic fields 0, 0, |], 
spontaneous nucleation of topological defects in expanding universes and false vacuum 
decay with time dependent initial states or time dependent potentials 0. In some special 
cases, these systems can be treated by standard instanton or imaginary time path methods; 
however, these techniques have limited applicability, and confusion often arises when one tries 
to extend the analysis to more general time dependent situations. For a discussion of various 
difficulties, see e.g. 0. 

In this paper we will investigate a collection of generalized time dependent versions of "stan- 
dard" tunneling problems, where the textbook instanton and imaginary time path methods are 
inapplicable due to the additional time dependence. The models typically have Lagrangians 
with an explicit time dependence arising from external backgrounds, or involve more compli- 
cated non-static initial states. Their unifying aspect is that they all can be analyzed via a 
method that combines the use of the WKB approximation with solutions of the classical equa- 
tions of motion along complex time paths. We will present a straightforward algorithm which 
can be used to compute the relevant tunneling or nucleation rates for such systems. 

To give a concrete example of this method, let us consider pair creation by a spatially 
constant electric field. In order to identify the specific features associated with a time dependent 
field, it is useful to first review the simple case of a static field. This problem, first solved by 
Schwinger [0, is most elegantly treated by an instanton approach. Calling the state with no 
particles present the false vacuum, the decay rate is determined by the imaginary part of the 
false vacuum energy, which can be extracted from an imaginary time path integral over fields 
which approach the false vacuum at r = —it = ±cx). In the semiclassical approximation 
one saturates the path integral by a solution to the (euclideanized) equations of motion; the 
corresponding configuration is the instanton. The action of the instanton determines the decay 
rate: T oc 6""^'"='=^"*°". 

The instanton solution describes the nucleation of a particle, anti-particle pair in the back- 
ground electric field. This can be seen directly by cutting the instanton in half. Half of the 
instanton solution corresponds to interpolating between the false vacuum at r = — oo, and a 
turning point, which we can take to occur at r = 0. The turning point configuration is that of a 
pair of particles momentarily at rest. If we were to continue evolving in imaginary time towards 
r = oo, then the particles would would converge and disappear, leaving the system in the false 
vacuum again. Instead, however, we can continue the solution to real time at the turning point, 
in which case the particles accelerate away from one another. Thus the full production process 
can conveniently be described by a combination of real and imaginary time evolution. 

For our purposes it is actually more convenient to consider the preceding discussion in the 
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reverse order. We can start by considering the real time expanding solution, and then consider 
evolving it back in time. Eventually we will reach the turning point, at which point we continue 
the evolution to imaginary time. If the particle separation proceeds to smoothly shrink to zero 
size in imaginary time, then the trajectory considered corresponds to a pair production process, 
and its action determines the decay rate. So, to summarize in a way that is most useful for the 
proceeding discussion, we look for expanding solutions, which can be smoothly shrunk to zero 
size when evolved back along some complex time contour. 

Phrased in this way, it is apparent how to adapt the procedure to the more general problem 
of a time dependent electric field. We can again look for solutions describing expanding pairs, 
but this time the continuation to complex time is more involved. Due to the time dependence 
in the problem, we no longer expect that the time contour along which the pair shrinks to zero 
size is one involving periods of purely real or purely imaginary time evolution; instead, the 
contour will be a more general curve in the complex time plane. Given that we can find such 
a contour, we can proceed to evaluate the action to determine the decay rate. The result will 
be a decay rate with non-trivial time dependence. 

By itself, the pair creation problem in a time dependent electric field has a long history. 
It was studied rigorously by Brezin and Itzykson 0, using Schwinger's proper time approach. 
Marinov and Popov treated it as a barrier penetration problem and employed WKB meth- 
ods. Further, as we discuss in the Appendix A, Audretsch noticed that the problem is 
isomorphic to overbarrier scattering in quantum mechanics. However, an advantage of the 
approach outlined above is that it is easily adapted to other problems involving the decay of 
metastable states via the production of extended objects. Further, it yields an instantaneous 
pair production rate with an explicit time dependence, so that one sees a time modulation in 
the flux of produced particles. The result for the pair production rate in applies only at 
specific times. 

As an example of adapting our approach to other physical processes, one can consider 
generalizing the computation of false vacuum decay in field theory, which proceeds through the 
nucleation of bubbles of true vacuum, to include field potentials with explicit time dependence. 
A particularly interesting source of time dependence arises from expanding universes, where it 
is expected that the expansion gives rise to the spontaneous nucleation of monopoles, strings, 
and domain walls. The nucleation rate has been computed for a very specific case, namely De 
Sitter space; but this is not in fact a time dependent problem, as the De Sitter geometry is 
static. For non-static geometries the more general approach discussed above is required. 

In the next section we write down an action which is general enough to treat the various 
processes we have referred to, and then give an algorithm by which one can compute the time 
dependent nucleation rate of the corresponding objects. In most cases, several steps in the 
procedure must be performed numerically. The simplest case, in which almost everything can 
be done analytically, is pair production in a time dependent electric field. We perform these 
steps in section 3, showing that the production rate takes a compact integral form. A special 
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case of this formula was derived before in 0. We review the connection of the problem to 
the problem of above barrier scattering in quantum mechanics. This discussed in detail in 
Appendix A. We then study the particular example of a sinusoidally varying field, and analyze 
the instantaneous pair production rates. In section 4 we turn to the other source of time 
dependence mentioned above, arising from the initial state rather than from external sources. 
Specifically, we consider a field theory with a local, but not global, minimum, and take the 
initial state to be one in which the field is undergoing coherent oscillations about the local 
minimum. We are able to calculate analytically for small oscillations, and to obtain the leading 
correction to the decay rate. Some computational details are relegated to Appendix B. Finally, 
in section 5 we summarize our conclusions and discuss directions for further study. 

2 Time Dependent Tunneling 

In this section we discuss our general approach to tunneling in models with explicit time de- 
pendence. Our aim is to show, using complex time contours, how such systems can be treated 
by a natural extension of the standard instanton method. To start with, we assume that there 
is an underlying field theory description of the model under consideration, and that that the 
system is initially in some m etastable state. We further assume that the state can decay via 
quantum tunneling, and that the decay occurs through the production of objects which can be 
described by a first quantized action. A wide variety of such objects can be described by the 
following action: 



Some specific cases of interest are: 

• a{x,t) = m ; b{x,t) = qE{t)x 

This is the relativistic action of a particle of mass m and charge q moving in a time dependent 
electric field E{t). As will be discussed in detail, this is the appropriate action for considering 
pair production due to the electric field. 

• a{x,t) = ATTa{t)x'^ ; b{x,t) = |7rp(t)a;^ 

This is the action of a spherical "bubble" of radius x, with time dependent surface tension 
a{t) and bulk energy density p{t). It describes false vacuum decay in a field theory from a 
time dependent initial state or in a time dependent potential, in instances where the thin wall 
approximation is valid. 

• a{x,t) = mc{t) ; b{x,t) = 

This is the action of a massive particle moving in the metric 

ds'' = c^{t)[dt^~dx''-s^{x){de^ + sin e^d(l)^)] 
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at fixed 9 , Here the choices s{x) = sinx, x, sinh a; give closed, fiat, and open Robertson- 
Walker universes. 



• a{x,t) = 2'7ific^(t)s{x) ; b{x,t) = 

We then have the action of a circular cosmic string moving in the above metric. The string is 
located at = 7r/2 and is centered at x = 0. 

• a{x,t) = A'Kac^{t)s'^{x) ; b{x,t) = 

This, similarly, describes a spherical domain wall in the above metric, again centered at x = 0. 

One might question the applicability of first quantized action to describe these systems, 
which are fundamentally field theories. In the electric field example, it is possible to make the 
connection rigorous and explicit, as is discussed in Appendix A. We can see no reason why the 
connection should not be valid in the other examples as well. 

We now explain how the action (|l|) can be used to compute the spontaneous creation rate 
of the objects which it describes. To begin, we should find the classical, real time, trajectories. 
The equations of motion are 



where ' denotes d/dx. For the purposes of tunneling, the relevant trajectories are those which 
emanate from a turning point. For a trajectory x(t), the existence of a turning point at time 
tf means that the canonical momentum vanishes there: 



In most cases, this condition will be simply a; = 0. A nucleation process corresponds to a 
trajectory which smoothly shrinks the object down to zero size when evolved back in time 
along a complex time contour. The continuation from real to complex time occurs at the 
turning point. Shrinking to zero size means that x = [| in the flat space examples, whereas in 
curved space it is the physical size c(t)s{x) which is required to go to zero. The condition that 
the shrinking to zero size be smooth is most easily seen in the flat space examples. Then we 
require that i; — ^ when x ^ 0; otherwise in the electric fleld case, for instance, the joining of 
the particle trajectory, x(t), and the anti-particle trajectory, —x{t), will be singular. In curved 
space a slightly more detailed analysis is necessary, depending on the speciflc form of c(t). To 
summarize, the trajectories of interest satisfy p(tf) = , x(to) = , x(to) ^ oo, tf = real, to = 
complex , x{t) = real. The problem is then: given some time tf, flnd an initial size x{tf) and 
complex time to such that the above conditions are satisfled. 

-"^In the electric field case, x denotes the position of a particle whose anti-particle is located at —x. x — 
thus corresponds to zero separation. 





(2) 




0. 
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Since x is required to be real while t is complex, the easiest way to proceed is to rewrite the 
equation of motion (H) as an equation for t{x): 

= -a'(x,t)Vt'2-l + 6'(x,t)t'. (3) 

The conditions on t{x) become 

= ; t'{0) = 0. 

*/ 

The advantage of this form is that it is straightforward to solve (|^) numerically, even if it is not 
possible to do so analytically. Then we can search for an initial coordinate value Xf = x{tf) 
which leads to t'{0) = 0. Having found an appropriate tunneling trajectory, t{x), we can proceed 
to evaluate its action. Again, it is easiest to change variables in (0), yielding 

S = J dx [— a(x, t{x))^Jt'^{x) — 1 + b{x, t{x))t'{x)]. 

The limits of integration run from to Xf simply because that gives the entire contribution to 
the imaginary part of the action. The imaginary part of the action determines the decay rate: 

Notice that the decay rate depends on the time tf which is the end point of the complex time 
path. Therefore, the result really is an instantaneous decay rate: the decay probability per 
unit time (per unit volume) at a time tf. Naturally, this is an idealized situation, a more 
realistic decay rate would involve a time average depending on the characteristic time scales 
and physical limitations of the given situation. However, even after such averaging, the decay 
rate characterizing the production rate will typically still be time modulated, and the time 
dependence can lead to observable consequences. 

3 Particle Creation in a Time Dependent Electric Field 

As a concrete illustration of the method discussed in the previous section, we will compute 
the probability for production of charged particles in a time dependent electric field. As noted 
previously, a rigorous approach to this problem involves starting from the second quantized 
field theory. For instance, analyzed the pair production in an alternating electric field 
through a Schwinger proper time approach. An alternative calculation could proceed through 
a Bogoliubov transformation relating the 'in' and 'out' asymptotic states; this approach is 
outlined in Appendix A. However, as we also discuss in Appendix A, the same results can be 
obtained from the more intuitive first quantized approach which we consider in this section. 



d 

dx 



a{x, t) 
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We begin by considering the simple case of a static field. That a (spatially and temporally) 
constant electric field should be capable of creating charged pairs is evident on energetic grounds, 
since the requisite energy, 2mc^, needed to create the pair is supplied by the electric field, if 
we separate the particles by a distance 2mc^/qE. The creation rate is proportional to the 
probability of the particles, initially located at the same point, to tunnel to this separation. 
Our starting point is the spin-0 point particle action in the presence of a constant electric field: 



S^-j dt[m\/l - i;2 - qEx] 



(4) 



For simplicity, we have taken space to be one dimensional. For ease of comparison with the 
discussion in Appendix A, it is actually simpler to integrate by parts and use the action 



S^- j dt[mVl - + qEtx] . 



(5) 



(This amounts to a different gauge choice for the gauge potential ^4^.) The equation of motion 
is 

d 



m— 
dt 



x 



x^ 



qE 



The solution is 



m 



x{t) = x{to) + — { 



\ 



1 + 



— {t-to) + 
m Ji 



i-2(to) 



Since we are interested in a tunneling process, we continue these trajectories to imaginary time 
t I— > —ir. Let us consider the particular trajectory for which Tq — —m/qE; x{to) — 0; x{to) — 
ioo. This yields the circle 

x^{t)+t^ 



m 



'qE' 



a quarter of it is shown in Figure 1. 



x=m / qE 



T=0 



x=0, T =-m / qE 

Figure 1: A trajectory for half of the tunneling process. 
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The particle comes to rest at r = after having traveled a distance Ax = m/qE. This 
portion of the path (one quarter of the circle) describes half of the tunneling process. The 
action for this portion is found by substituting the trajectory x(r) into 



-ra/qE 

where ' now means d/dr. This yields 



S = i C dr[m\ 1 + (-^)2 + qETx(T)] 

J^mlqE V dr 



S 



AqE 

The full tunneling process is described by a semicircle (see Figure 2.) 

X = -m / qE X = m / qE 



T = 



X = , x= -m/qE 

Figure 2: A trajectory for the full tunneling process. 

Whereas the part of the trajectory between x = and x = m/qE depicts the creation of 
a particle with charge g, the part between x = and x = —m/qE depicts the creation of the 
corresponding antiparticle of charge —q. We see that the antiparticle trajectory is obtained 
from the particle trajectory by g —q; x(r) — x(r). It is easy to see that this path yields 
a solution for the equation of motion and has the same action as the particle trajectory. The 
total tunneling action is thus 

inm"^ 
2qE 

The rate of pair creation is found by squaring the tunneling amplitude: 



g-Im[Stotal] 



2 7rm 

= e 



The exact result obtained by Schwinger is 

\2 oo / 1 \n+l 



{qEf 
(2vr)3 

The WKB result gives a good approximation when irm'^/qE ^ 1. Finally, we note that after 
being created, the particles move along hyperbolas, as is seen by continuing the trajectories 
back to real time. 

We now generalize the analysis to treat a time dependent electric field. As discussed in 
the previous section, we expect that this can be accomplished by finding a complex time path 
connecting the initial and final positions. 
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In particular, we will look for a path as shown in Figure 3. 




X = 0,t= 



Figure 3: A trajectory for half of the generic tunneling process. 

which describes the creation of a particle at time t — tf. The vertical axis no longer represents 
purely real or imaginary time, but rather some more general complex time direction (such that 
tf is real). At t — tf the particle is at rest, x{tf) — 0. At t — tq the velocity should be 
singular, x{to) — > oo, so that the trajectory for the antiparticle can be smoothly joined on to 
the particle trajectory. Both of these conditions were of course satisfied in the static field case. 
The remarkable aspect of this problem, as we shall see, is that these conditions allow us to 
determine the tunneling action without having to find the complex time path explicitly. 
The action, after integrating by parts, has the form 

S^- J dt[mVl - - qA(t)x] . (6) 

where we have defined ^ 

A{t) = - I dt'E{t'). (7) 
Jtt 



The equations of motion then yield 

x{t)^ 



qA{t)/m 



l + q^A^{t)/m? 



We have set the conserved momentum — dL/dx to zero, so that x respects the condition 
x{tf) — 0. Integrating: 



x(t) = /* dt' , ^^^'^ + xitA 

^ Jif . l + Q^A^{t')lw? 



+ q^A^{t')/w? 

Now, to is determined by requiring i;(to) oo. For a nonsingular E{t), this implies A{tQ) 
±im/q. We choose to set A{tQ) — —im/q. 
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We will now write the action in terms of A{t) and see that it takes a simple form. Substi- 
tuting the expression for x into (Bf) gives 



S = — m 



dtJl + q^A\t)/m^ 

Jto ^ 



Now we change the integration variable from t to v = —qA/m using 

m dv 



dt 



qE{t{v)) 

where t{v) is found by inverting (0). Then, 



S = / (if—-—- = / dO- 



q Ji E{t{v),tf) q Jo E{t{isme),tf) 

This is the action corresponding to the creation of the charge q particle. As before, the action 
for the antiparticle is obtained by replacing q —q; x{t) —x{t), which yields the same 
result as for the particle. Therefore, the pair creation rate is given by 

r(t;)=exp{-41m[5(t^)]} . 

To recapitulate, the essential trick that was used was to map the potentially complicated 
complex time contour into the complex f-plane, where it always takes the simple form of a 
line from to i. This simplifies the problem considerably, since it is no longer necessary to try 
to find the complex time contour. We only need to know its image in the f-plane, and we do. 

3.1 Example of Time Dependent Pair Creation 

Having obtained the general result^ 

im'^c^ cos^6' , , 

S = / dO- 



q Jo E{t{icsm6),tf) 
for the action along the complex time path, we shall now evaluate the pair creation rate 

T{tf)=exp{-^lm[S{tj)]} 

in an example case of a time dependent electric field. But first, we check the formula with a 
constant electric field E = Eq. Now the denominator in (|^) is a constant, and what remains is 
an elementary integral, so we easily obtain the standard WKB result for the pair creation rate 

2 S 



^In this subsection, we restore U and c. 
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This pair creation rate is generally truly small. Even for strongest electric fields obtained in 
a laboratory, Eq ~ 10^^ N/C |§, the exponent is still enormous: using g = e ~ 10~^^ C, 
mc^ ~ nieC^ ~ 10^ eV ~ 10^^'^ Nm, c ~ 10^ m/s, h ~ 10~^^ Js, we get a vanishingly small 
production rate 

r ~ exp{-10*^}. 

Now, let us consider an example of a time dependent electric field, an oscillating strong 
electric field 

E(tf) = Eq cos tut/ . 

Note that no magnetic field is required to satisfy the Maxwell's equations, only an alternating 
uniform current density. This case is different from pair creation in the background of an 
electromagnetic plane wave which was considered originally by Schwinger in the sense that 
the intensity of the plane wave is time independent. 
We first find v and E{t{v),tf): 



v(t,tf) = — [ dt' Eit') = - — -( sin ut — sin utf 
m Jtf mio 



and 



^/ / N ^ / .iraucsinO 
E{t{icsme),tf) = Eo\ I - ( + smwt/)2 

V 



After changing the integration variable to a; = sin 6, the action can be written as 



,2„3 ^1 



ij — — zz — / (lX~ 



qEo Jo ^1- (i^a; + sinwt/)2 

Finally, isolating the imaginary part of the action we find the instantaneous pair creation rate 
to be of the form 

where /(t/) is a modulation factor which characterizes the time dependence (compare with (^). 
We find that 

2^2 /■! Vl-x^Jr{x,tf) +uix,tf) 

litf) = / dX r , 

n Jo r{x,tf) 



u{x,tf) = 1 + x'^ — sin^ tot f 



r{x,tf) = ^Ju'^{x,tf) + Aa'^x'^ sin ujtf 
mujc 



a = 



qEo 

This integral can be evaluated numerically, we show a plot of /(t/) for a = 1 (u; = qEo/mc): 
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1 ' ' ' ' i ' ' ' ' i ' ' ' ' i ' ' ' ' t ' ' ' ' i 



Figure 4: An example plot of the time modulation function I{tf) in the action. 

We can see that I{tf) is a periodic function of time. It reaches a minimum value at times 
utf = + mr and a maximum at utf = 7r/2 + nvr. The dependence of the minimum and 
maximum values of / on the frequency lu turns out to be interesting: 




Figure 5: A plot of the minimum and maximum values of the action as a function of 
frequency. 



At zero frequency, the minimum value of the time modulation function is equal to 1. The 
maximum value becomes infinite: at small frequencies the decay rate is completely dominated 
by the minimum value of the action. Thus at zero frequency the decay rate reduces to that 
of the static it should. As the oscillation frequency increases, the maximum and 

minimum values decrease monotonically and both seem to approach the asymptotic value zero. 
This means that the (average) pair creation rate increases as the oscillation frequency of the 
field increases. Finally, the action becomes so small (and the rate so high) that the WKB 
approximation is no longer valid. 

Let us try to check the asymptotic behavior of the maximum and minimum. For the 
minimum, this can be done rigorously. Setting sin cut/ = simplifies the integral and we can 
identify it in terms of complete elliptic integrals. We find 

41/1 ,f^^/ CL , ^, a ,1 
Ti y ^a? + 1 Va + 1 J 

This result was also found in [^. As a; — >• oo (a — >• oo) we get 

In a 

Jmin ~ ^ 

a 

for the leading asymptotic behavior. For the maximum value, the analysis is a bit more com- 
plicated, but we find the same leading asymptotic behavior 

In a 

/max ~ ^ . 

a 
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4 Decay of a False Vacuum 



We now turn to another problem involving tunneling with non-trivial time dependence, whose 
treatment requires methods slightly different from those we have discussed to this point. Let 
us consider a scalar field 



Figure 6: Potential. 

We take the initial state to be one in which the field is concentrated in the well centered 
at = 0/. The field will be taken to be constant in space, but can have a non-trivial time 
dependence. In particular, we have in mind semi-classical looking states in which the field 
oscillates coherently in the well. Such a configuration, though stable classically provided the 
amplitude of oscillation is below the barrier, is expected to be unstable quantum mechanically. 
The presence of the well at (p = (pt signals a decay process whereby the field can tunnel through 
the barrier. This process is described by bubble nucleation, meaning that regions of field 
concentrated at (pt spontaneously form within the initial configuration and rapidly expand. 
We would like to know the rate at which bubble nucleation occurs, and describe the resulting 
bubble trajectory. It is also important to determine the state of the field inside the bubble - we 
might imagine that that the oscillations about the false vacuum outside the bubble feed into 
the interior of the bubble, causing the field there to oscillate about the true vacuum. By solving 
the field equations we will see that such oscillations are actually confined to a region near the 
bubble wall, so that the interior field in the bulk of the bubble is frozen at the true vacuum. 

We begin by considering the simplest case, where the field is initially located at the bottom 
of the leftmost well, (f){t) = (pf. This, of course, is the case considered by Coleman P]. Our 
strategy will be to look for an expanding bubble solution which can be shrunk to zero size when 
evolved back along a complex time contour. This problem is most efficiently solved by utilizing 
the SO (3,1) symmetry of the theory. However, we will later be considering initial states which 
oscillate coherently, and break the S0(3,l) symmetry to S0(3). Therefore we will discuss the 
solution in a language which explicitly uses only the latter symmetry. 




(10) 



and take the potential to be of the form shown in Figure 6. 
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For simplicity, wc will work in the thin wall approximation, which is valid provided the 
difference in energies between the true and false vacua, p = — V{(f)t)-, is sufficiently small. 

The bubble solution then has the form of a spherical region of true vacuum separated by a thin 
wall from the outside region of false vacuum, 

<t>Kr,t)^^ 0^ for r>i?(t) . ^^^> 

The trajectory of the bubble wall, R{t) , can be determined from energy conservation. Consider 
the energy in the region r < R. There arc two contributions to the energy. The interior of the 
bubble contributes £'(inside) = ^nV{(f)t)R^ ■ There is also an energy proportional to the area of 
the bubble wall associated with the field gradient in passing from true vacuum to false vacuum: 



£'(wall) — 47r(7o-R^/Y 1 ~ Here, uo is the energy density of the wall. 



(To = / dr 

J wail 



in the next section we will compute its value from the field equations. The energy £J(inside) + 
£^(wall) must be equal to the energy present in the region before the nucleation of the bubble: 
£^totai = lnV{(j>f)R^ So, 

47rc7oi?2 4 



- -Trpoi? = (12) 



with 

The trajectory is then 
with 



p,^V{<f>f)-V{<f>t). 



R{t) = ^JRl + 



i?o = ^ . (13) 
Po 

For t > this describes an expanding bubble solution. To consider tunneling, we evolve the 
solution back to the turning point at t = 0, and then try to shrink the bubble to zero size along 
a complex time contour. In the present case this step is trivial - the contour displayed in Figure 
7. does the job. 
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Im(l) 

T V"*0 



Reft) 



Figure 7: Time path for bubble nucleation. 

It remains to determine the amplitude for the tunneling process, and for this we require the 
classical action. The bubble has a Lagrangian 



Since we wish to compute the relative probability of bubble nucleation versus remaining in the 
false vacuum, what we actually want is the difference in action between the bubble solution 
and the false vacuum state. This is given by 

4 



S = Jdt Lbubbie + -7rl^(0/)^' =-/ 



dt 



3 



The action can be put in a useful form by inserting the equation of motion (12), and changing 
variables to R: 



The action has an imaginary part coming from the part of the trajectory < i? < Ro, when 



S = J dRR'^{^porR^-{Anaoy = J 

las an imag 
the bubble is tunneling: 



The nucleation rate is then 

r^e-^^-[^]=exp|-^Poi?o'j 

which is Coleman's result. 

Now let us generalize to the case where the field is initially oscillating around the false 
vacuum: (p = (j)f{t). In we will study the bubble solutions with this initial condition, and we 
summarize the results here. For small oscillations (j)f{t) = (pf + a{t) the bubble looks like 

■pt for r < R — A 

0bub(r, t)^{ 0t + f exp[(r - R)/A]a{t) hi R - A < r < R (14) 
(pfit) ioY r > R 
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where A is small compared to the typical size of the bubble R. In other words, the field 
oscillations only penetrate a relatively small distance into the bubble; the bulk of the bubble's 
interior simply sits at the true vacuum as before. As we have already mentioned, although one 
might have expected the bubble to leave a state oscillating about the true vacuum, we see that 
this is not the case. 

Now we invoke the same energy considerations as before. Since the field oscillations inside 
the bubble are localized near the wall, the energy inside the bubble is essentially given by the 
true vacuum configuration 0^. Thus 

4 

-E^bubbie (inside) = -W {(t)t)R^ ■ 

The bubble wall has the energy 



-Sbubbieiwall) = 



where 



dr U(<^bub)' + ^(0Lub)' + n0bub)} . 



/ wall L 2 2 

(Te is time independent. The initial energy is also divided into two contributions: 



-^initial (inside) = --kR^ 



and 



where 



(wall) 



initial 



'E - dr 

^wall 



FV 



Conservation of energy then requires 

AiiaER^ 



R^ 



i?2 



.{4,f{t)f + v{4>f{t)) 



- -ttpeR^ 



with 



O-E 

Pe 



-bubble 



FV 



'E 



(15) 
(16) 



We emphasize that pe and as are constants. This fact means that the complex time contour 
relevant for tunneling runs in the purely imaginary direction, just as in the static case. However, 
we expect that this behavior is an accident of the analysis in the limit of small oscillations; 
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more generally, pE and o"e will acquire time dependence and the time contour will be a more 
complicated curve in the complex time plane. Now we define 

3(7e 



Rq 



Pe 



so that the trajectory is 

Kit) = ^Rl + {t- . (17) 

Now that we have obtained the bubble trajectory, we turn to the evaluation of the decay 
rate. As before, the decay rate is found by integrating the the action over an imaginary time 
contour running from some initial time Iq, to + iRq where the bubble shrinks to zero size. 
The action which is to be integrated is the difference between the bubble action and the false 
vacuum action. The bubble action is 



^1 



' wall 



bubble 

where 

^bubble _ _ f 

J w 

The false vacuum action is 

Spy 

where 



- / dt 



ATiaY^'°^\t)R\t)^l - i?2 + ^V{(t)t)R 



47r, 



dr \-{ 



2(<^bub)^ 



\ 1 ^bubble f 
bubj f = CTfi - / 

) J wall 



dr (p] 



bub 



dt 


47rarR\ 


li- 




= 4^- 


1 

J wall 


pr 


= p7- 


0?. 



dr 



The action to be integrated is thus 



S = - J dt AT^aL{t)R^{t)\ll- - YPdt)R 



(18) 



where 



PL{t) 



dr 



wall 
2 



^bub 



(19) 
(20) 
their 



(^E 

Pe - 0/ 

A crucial point is that although as and pe are constants, ai and pl are time dependent 
time dependence determines the time dependence of the decay rate. 

The calculation has thus been reduced down to performing the integrals for a^, pl, and S. 
These are straightforward to do; they can be done analytically for small oscillations about the 
false vacuum, as is shown in Section 4.2, although in the general case numerical integration is 
required. The result is an expression for the time dependent decay rate: 



r(io) !^exp{-2Im[5(to)]} 
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4.1 Structure of the Oscillating Bubble 

To analyze the case when the field is initially oscillating around the false vacuum, (f) — (f)f{t), 
we first need to determine how the structure of the bubble is altered. Wc will now present 
an example calculation of the bubble solution 0bub('")^) which interpolates between the true 
vacuum (f)t and the oscillating initial state (/)f{t). We shall refer to this as the oscillating bubble. 
The field equation that we need to study is 

As an example, we consider the potential discussed by Coleman, 
where e > 0. The true vacuum is located at 

where the potential has the value ^(0-) ~ — e + 0{e^). The false vacuum is located at 

where ^(0+) ^O + Oie^). 

In the standard scenario of decay from false vacuum to true vacuum, the structure of the 
bubble is obtained from the static solution, (poit^), of the field equation, 

00 + ^00 = ^(00). (21) 

00 interpolates between the true vacuum 0_ inside the bubble, and the false vacuum 0_|_ outside 
the bubble, with the non-trivial r-dependence concentrated in the bubble wall. Without going 
into the mathematics of the exact form of the solution, we recall that the qualitative behavior 
of the solution is as depicted in figure 8. 
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Figure 8: Picture of the qualitative behavior of the static solution (l)o{r). 



The radius of the bubble is R = 3cro/e (see eqn. (pTSD) and the thickness L of the bubble 
wall is of the order L ~ where 



An approximation to the static solution in the vicinity of the wall can be obtained by dropping 
the term (2/r)0Q in the static field equation (since r ~ ^ 0), and dropping the small constant 
term e/2a. The solution of the approximate field equation is the 'kink' 

where yU = ay/X. Its behavior is similar to that depicted in figure 8. 

Now we try to find a time-dependent solution 4>hnh{r,t), which reduces to the coherently 
oscillating field 0/(t) = 0+ + sin tut about the false vacuum as r — > oo. We assume that the 
amplitude of the oscillations is small. The frequency u of the oscillations is given by 

.^.0(,,).4A„^_il,O(a. (22) 
We make the following ansatz for (f)huh{r,t): 

4>hnh{r,t) = 4>o{r) + a{r) sinut . 

Substituting this ansatz into the full field equation, and using the fact that 0o is the static 
solution, we obtain a linearized differential equation for the profile function a{r): 



2 

a"(r) H — a'(r) + 
r 



' 2 d-'V ■ 



a{r) = . (23) 



Without knowing the exact form of the static solution </)o, we can still proceed by using its 
known asymptotic properties. In the region outside of the bubble, r > R, (j)^ reduces to 0+, so 
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{cP'V / d(j?){(j)Q) — > {d^V/d(j)^){(j)^). We see that since the frequency u"^ is equal to 

the profile function a{r) must be constant in this region. So (phubif^i) correctly reduces to the 

oscillating initial configuration in this region: 

0bub(^, ^) ~* 0+ + «o sin ut as r > R . 

In the region inside the bubble, r < R, reduces to the true vacuum value 0_. Thus 
^ (rf2v^/rf02^(0_). Using 

and denoting A;^ = 3e/a^, we see that the differential equation (^) reduces to a familiar equation 

2 

a"(r) H — a'(r) — A;^a(r) = . 
r 

The solution for a{r) is, 

sinh(A;r) 



a(r) = A- 



kr 

We fix the constant A by matching a{r) with ao at r = i?: 

kR 

A = aQ — 



smh.{kR) 

The solution for a{r) tells us that the oscillations a{r) sincijt decay to zero inside the bubble, 
in a region of thickness A = 1/k. Thus there are three scales that characterize the structure of 
the oscillating bubble: 

1. The radius of the bubble R ~ ^co/e. 

2. The thickness of the bubble wall L ~ l/{a^/X). 

3. The thickness of the region inside the bubble where the oscillations decay A ~ a/{3^/e). 

The relative sizes of these scales are as follows: 

A ay/e 



R o-Q 

L 



A a^yX 

Thus, 



as e — s> 
as e — > . 



Also, 



L < A < i? . 



fci? ~ ^ » 1 
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Hence the solution for a{r) can be rewritten as 

r 

to a good approximation. 

Finally, we would like to obtain at least an approximate solution for a{r) in the bubble wall 
region r R. We do this by replacing 0o(^) in the equation with the approximate solution 
a tanh /i(r — i?), and dropping the term (2/r)a'(r). Further, we approximate the frequency 
by uj"^ = 4Aa^. Then the equation (^Sf ) reduces to 

6 

a"{x) H 2 — '^(^) = ) 

cosh X 

where x = n{r — R). This has the solution 

a{r) = — [3tanh^(/i(r - R)) - 1] . 

Since a(r) — 5 in regions r < R where we know that a{r) = a^, we set B = Oq- 

To summarize, we have found an approximate solution for the oscillating bubble: ^bubl'"; t) = 
0o(r) + a(r) sin cut, where (poii^) is the solution in the static case, modified by small oscillations 
with a profile function a{r) given by 

f ao r^R + L/2 

a{r) = I f [3 tsinh\{r - R)/L) -1] R- L/2 ^ r ^ R + L/2 (24) 

[a^Eeir-R)/A r^R-L/2. 

4.2 Calculation of the Bubble Nucleation Rate 

Now that we have found the oscillating bubble solution (f)buh{f,t), we can proceed to calculate 
the instantaneous bubble nucleation rate r(to) ~ {2Im[S'(to)]}, where to is the time immediately 
after the bubble has nucleated. 

First we need to find the quantities and pe that appear in the equation of motion of the 
bubble. They are: 

(To - 0.42 co^alL (25) 

,-'-2 2 

e + -uj ao . 

For some calculational details, see Appendix B. Since (Te,Pe are constant, we could use energy 
conservation which gave us the trajectory of the bubble (|T7|): 

R=^Rl + {t- toy , 



aE ~ 
Pe ~ 
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where Rq = Sas/pE- 

For bubble nucleation, we need a time path which shrinks the bubble to zero radius. Thus 
the time path has an imaginary segment 



t = to + ^^JRl- 

for R < Rq. The above branch choice gives the correct exponential behavior for the nucleation 
rate. 



To calculate the imaginary part of the action, it is first convenient to divide the action (18) 
into two parts: S = Sq + Si, where 

^0 = - r dt (47rasi?Vl-i?2 _ ^peR^\ 

Jtn+iRo I 6 J 



and 



rto \ I — /•R+L/2 , . . X Att . 

Si= dt\ 4vri?2vT^ / dr - - 

Jtn+iRo JR-L/2 ^ ^ 6 



ito+iRo y JR-L/2 

For the first term S'o, the calculation proceeds as in the 'static' case. The result is 



^2 



ImSo = —PeRq ■ 
For the second term, after a bit longer calculation we find 

ImSi ^ aln^Rl | ^^"^^^ + [0.84 h{2ujRo) + ^ l2{2ujRo)] cos(2c<;to)| • 

[We have used L ^ 2/u and dropped a negligible subleading term.] The functions /„(2u;_Ro) 
are modified Bessel functions. The total instantaneous bubble nucleation rate is then 



2 ^(2ujR )^ / 1 

^ °^ ' '1.68 /i(2cji?o) + - /2(2cui?o) ) cos(2cuto 



T{to) ^ exp <j -yPiji?o - a^Tc^Rl 
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The decay rate is oscillatory, with the leading correction to the static decay rate coming from 
the oscillatory term in the exponent. Recall that Rquj ~ radius of the bubble / thickness of the 
bubble wall, so Rqu ^ 1. Thus, the leading order correction is of the order a'^Rle^'^^ /2uRq: 
this is much larger than what one might have anticipated. However, this is reminiscent of what 
happens when a particle tunnels through an oscillating barrier. Biittiker and Landauer showed 
rU| that the tunneling particle absorbs quanta from the oscillating barrier; the net effect is 



that tunneling becomes easier. In the leading order correction to the tunneling probability, 
the amplitude of the oscillations of the barrier is multiplied by an exponential term, so the 
correction is much larger. As we can see in our case, small oscillations about the false vacuum 
also render the state more unstable. 
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5 Conclusion 



In this paper we have presented a rather general approach to treating time dependent tunnehng 
problems, and have illustrated the method with some concrete examples. Through the use of 
the WKB approximation, we found that we could reduce such problems to solving classical 
equations of motion along complex time contours. Even when such equations of motion are 
analytically intractable, we have shown that they are amenable to straightforward numerical 
analysis. The standard approach to tunneling problems in field theory, using instantons, is seen 
as a special case of the complex time method. In particular, the familiar procedure of evolving 
along the imaginary time direction is valid when there is no non-trivial time dependence in the 
problem, but more general problems require more general complex time contours. 

The most straightforward and elegant application of our methods is to the case of pair 
production by a time dependent electric field. As discussed previously, this problem has been 
analyzed by a number of workers over the years by a variety of methods; we derived a slightly 
more general result than had been previously been obtained, but our main interest in the 
problem was as a prototype for more complicated time dependent systems. The qualitative 
features of the electric field problem carry over to these systems, the only difference being that 
several steps must be performed numerically rather than analytically. 

As an example of an interesting process which can be tackled by the methods developed here, 
we mention again the quantum nucleation of defects in an expanding universe. The form of the 
expansion can have an important effect on the resulting distribution of the defects. Consider 
cosmic string nucleation, for instance. The strings nucleate with a size equal to the horizon, 
or inverse Hubble constant, which thus also sets the scale for the duration of the complex time 
evolution. If the expansion rate of the universe varies appreciably on this time scale, then the 
nucleation rate will depart from what a naive quasi-static analysis would indicate. Depending 
on the cosmo logical model, the resulting distribution of strings is potentially relevant. 

To illustrate another source of time dependence - that arising from initial conditions, rather 
than external sources - we considered the problem of false vacuum decay in field theory, where 
the initial state consists of coherent field oscillations about the false vacuum. Again, instanton 
techniques are not directly applicable to this system. Although the analysis was rather involved, 
we were able to obtain an expression for the time dependent decay rate in the case of small 
oscillations. Of course, to see large time dependent effects one must allow for large oscillations, 
but this would require a rather intricate computation which we have not attempted. But again, 
we stress that the steps are, in principle, straightforward. 

Appendix A 

In this Appendix, based on ref. 0], we discuss how the problem of pair production in a 
homogeneous time dependent electric field is related to overbarrier scattering in nonrelativistic 
quantum mechanics; this connection is the basis for the first quantized approach presented in 
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Section 3. 



We start from a relativistic field theory formulation, in which a charged scalar in an electric 
field has the equation of motion 

{{d^ + iqA^){d^' + iqA>') +m'^} (f) = . (26) 

In the canonical quantization approach, we seek 'natural' mode solutions of (p6|) with suitable 
asymptotic properties, which we then use in the oscillator expansion of the field operator to 
identify the asymptotic states. If we consider the background of a classical time dependent 
homogeneous electric field E = E{t)ex, it is convenient to take the gauge choice A^^ = {A^, A) = 
(0, A{t), 0, 0), where A{t) = — / dtE{t). Then we can use the separable ansatz 

in the field equation, to reduce it to a time dependent generalized oscillator equation 

/p^ + /p^=0 (27) 

where uj'i{t) = rn? + {p^ — ^A{t)Y + + = /^^ + {Vx — eA{t)Y. In addition, the relativistic 
Klein-Gordon scalar product reduces to a simple form 

(/l,/2)=^(/l72-A72) . (28) 

Equation (^) is identical to a nonrelativistic time independent Schrodinger equation 

V'" + A;2(x,E)V^ = , 

via the following identifications: 

t^x , f^^^lj (29) 
ul{t) =i? + (p.^ - qA{t)f ^ A;2(x, E) = 2mE - 2mV{x) . (30) 

Also, the scalar product (|28|) can be seen to match with the inner product for wavefunctions 
ipi,ip2- If we make the more precise identifications 

/i^ ^ 2mE (31) 
{px-qA{t)f ^ -2mV{x) (32) 

we notice that in the quantum mechanics language we must take E > and < 0. This 
suggests that we are dealing with a quantum mechanical overbarrier scattering problem. 

To complete this connection, we first outline the pair production calculation in the field 
theory picture. If we switch the electric field on in the far past and off in the far future, the 
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gauge potential is asymptotically constant: A{t) ^m.out ^ _^ ±00. Then the oscillator 
equation has two natural linearly independent solutions with the asymptotic properties 

t), t ^ ±00. These give two bases for the oscillator expansion of the 



exp( 



in, out 
P 



field operator, corresponding to the 'in' and 'out' Fock spaces, as in a quantum field theory in 
curved space [0. In particular, the definition of an initial vacuum state is different from a 
final vacuum state; this gives rise to particle production. To compute the particle production 
rate, one must find a Bogoliubov transformation which relates the two bases: 



jp 



rout _i_ n n*out 
^pJp ~r PpJp 



(33) 



Then the time averaged pair production probability F is given by 

2 



Pp 



dp 

To see the connection to the quantum mechanical scattering problem, we first rewrite 

1 

Oip 

Then, asymptotically as t ^ —00, the l.h.s. reduces to 



3f) as 



out ^ 

Op* 



*out 
P 



Jn 



-ioj^^'t I f^p iujZ^ 
g p _| g p 



ikx 



+ Re 



ikx 



Cif. 



and the r.h.s. reduces to (as t 00) 



Te 



-ik X 



We see that the analysis of the oscillator equation (^) in the field theory picture exactly 
corresponds to an overbarrier scattering problem in a quantum mechanical picture, with the 
refiection and transmission coefficients being related to the Bogoliubov coefficients by 

Pp 



\R\ 



\T\ 



The time averaged pair production probability is 

T = \R? 



If the electric field E{t) is not very rapidly varying, we can simplify the problem further with 
the use of WKB approximation. We use the WKB solution of the oscillator (or Schrodinger) 
equation 

/.expH/.*}"exp{-» (34) 
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The treatment of semiclassical reflection above a barrier requires an extension of the WKB 
One method is to first find the complex turning points off the real axis where 
0. Then, the leading contribution to the reflection coefficient can be found by 



method |T2 , 
a; ~ /c 

computing the integral / ujdt along a complex time contour C, traveling from some initial time 
tl on the real axis to the closest complex turning point to and back (see Figure 9.). 



,Im(t) 



Re(t) 



Figure 9: The integration contour C. 
This procedure yields the reflection coefficient 

= exp{-21m I to dt} . (35) 



Alternatively, we can derive the exponent of the WKB wavefunction ( P^ by starting from the 
action 

S = -j dt {mVl - x2 _ qA{t)x} (36) 

for the relativistic charged particle in the background electric field (using the same gauge choice 
as before). Using the equation of motion 

X qA 



and substituting x, gives 

S = -J dtVlT^ 

= - J dt^J{px - qAy + w? = - J dt u . 

The WKB wavefunction ( |5^ ) is e*'^ as expected. Therefore, we can use the one-particle action 
to calculate the reflection coefficient (^); this is the basis for the first quantized analysis of the 
pair production rate presented in Section 3. 
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Appendix B 



We now calculate the quantities and pe that appear in the equation of motion of the 
bubble. The energy density pe is: 



PE = ^(0/(t)f + 



-uj'a'o cos'ujt + V{(f)+) + - 



l<fV 



2 2 I 

ag + e 



On sin^ ojt + e 



For the surface tension, we first need 



bubble 



dr 

wall 1 2 



1 , • 



bub)^ + -0bub + ^(0bub) 



Substituting (phnhij^t) = (po{r) + S{r,t) = 4>o{r) + a (r) sin cut and expanding to second order 
in 6, we find that the leading order term gives the surface tension (Tq in the traditional false 
vacuum decay, 

dr 



wall 



The contribution from the linear order in S vanishes after integration by parts and using the 
equation of motion for 0o- The contribution from the second order in 6 is 



= - 



2 J wall 



dr 



'd^V 



a 



^>o) -uj']+ [a 



i\2 



• 2 . I 2 2 

sm uji + uj a 



Using the differential equation (23) for a (dropping the 2a' /r term) and integrating by parts, 
we find 

^2 i.R+L/2 



0"2 



2 Jr~L/2 



dr a (r) 



Substituting a{r) = ^[3tanh^((r — R)/L) — 1], and combining the contributions. 



bubble ^ 

U TP ~ 



ctq + 0.16 -io'^alL . 



There is also the contribution from the initial state 

1 



a 



FY 



dr 



wall 



^ 2 2 r 

2 ° 



Then, finally, the surface tension is 



bubble ^FV ^ ^ n /lO , ,2^2 r 

— cr^ ^ ctq — U.4/ uo a^L . 
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